function thebij(term,type)
load hermitdata;

if(strcmp(type,'polychaos'))
for jj=1:term
    for kk=1:term
        bij{jj}{kk}= @(x) hermitijk(1,jj,kk)+hermitijk(2,jj,kk)*x/sig;
    end
end

elseif(strcmp(type,'naive'))
bij{1}{1}=@(x) x;   bij{2}{1}=@(x) x;  bij{3}{1}=@(x) 2*x;   bij{4}{1}=@(x) 3*x; 

bij{1}{2}=@(x) 2*x; bij{2}{2}=@(x) x;  bij{3}{2}=@(x) x;     bij{4}{2}=@(x) 3*x;

bij{1}{3}=@(x) 3*x; bij{2}{3}=@(x) x;  bij{3}{3}=@(x) 3*x;   bij{4}{3}=@(x) 3*x;

bij{1}{4}=@(x) 2*x; bij{2}{4}=@(x) x;  bij{3}{4}=@(x) 2*x;   bij{4}{4}=@(x) 2*x;
end
fname='funcbij';
save(fname,'bij');